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By employing a numerical renormalization group technique, we evaluate electronic specific heat coeffi- 
cient 7 of the Anderson model coupled with local anharmonic phonons for the oscillation of a caged atom. 
For the rattling-type cage potential with a flat and wide region in the bottom, we find that phonon-mediated 
attraction is largely enhanced. When the potential shape is deformed from the rattling type, there occurs 
a cancellation between Coulomb repulsion and the phonon-mediated attraction. In such a situation, spin 
and charge fluctuations are comparable to each other, leading to the realization of exotic electron-phonon 
complex state with large and magnetically robust 7. 
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Recently, a new possibility of electron-phonon coupled 
state emerging in strongly correlated electron materials such 
as clathrate compounds, 1 pyrochlore oxides, 2,3 and filled 
skutterudites 4 has attracted much attention in the research 
area of condensed matter physics. A common feature of these 
materials is the existence of nano-size cage composed of rela- 
tively light atoms, in which another atom feels a highly anhar- 
monic potential and oscillates with large amplitude. Such an 
oscillation is frequently called "rattling", which is considered 
to be one of key ingredients of cage-structure materials, when 
we attempt to clarify their electronic properties. 

For instance, in order to understand magnetically robust 
heavy-fermion behavior observed in Sm-based filled skutteru- 
dite compound SmOs4Sbi2, 5 the non-magnetic Kondo effect 
originating from phonon degree of freedom has been pointed 
out. 6 8 Along with this research direction, in order to promote 
our understandings on the Kondo physics in electron-phonon 
systems, the present author has performed numerical calcula- 
tions of the Anderson model coupled with local Jahn-Teller 
and Holstein phonons. 911 The periodic Anderson-Holstein 
model has been also analyzed with the use of a dynamical 
mean-field theory and a mechanism of the mass enhancement 
due to large lattice fluctuations and phonon softening towards 
double- well potential has been addressed. 12 This tendency has 
been also suggested in the Holstein model. 1315 

Concerning the robustness of electronic specific heat co- 
efficient 7 against an applied magnetic field, it is basically 
understood, if the Kondo temperature Tj< becomes large in 
comparison with the magnetic field. In fact, in an electron- 
phonon system, Coulomb interaction is considered to be, 
more or less, weakened by phonon-mediated attraction, lead- 
ing to the enhancement of 7k- However, in order to cancel 
the on-site Coulomb interaction, phonon energy should be in 
the order of electron bandwidth, as long as we consider har- 
monic phonons. Such a situation may be realized in some 
molecular conductors composed of light atoms, but in gen- 
eral, phonon energy is smaller than the electron bandwidth. 
Moreover, when we recall a relation of 7 ~ 1/Tk, 16 7 cannot 
be so large for the case with enhanced Ta- 
in this paper, in order to examine the phonon-based sce- 
nario for large and magnetically robust 7, we investigate the 
effect of anharmonic oscillation of the atom in the cage poten- 
tial by treating them as local phonons in the standard quan- 



tum mechanics. We analyze the Anderson model coupled 
with local anharmonic phonons with the use of a numerical 
renormalization group method. For a rattling-type potential 
in which a flat and wide region appears in the bottom, the 
Coulomb repulsion is largely suppressed by phonon-mediated 
attraction, even when we consider phonon energy smaller than 
the electron bandwidth. We also find that 7 is both large and 
magnetically robust, when the potential shape is deformed 
from the rattling type in which spin fluctuations are compa- 
rable to charge ones. 

First let us discuss the local electron-phonon state, when 
we include anharmonicity. The local Hamiltonian is given by 

H\ oc = Un^n l + H cph , (1) 

where n CT =dJ.d CT , d a is an annihilation operator of localized 
electron with spin a, U is the Coulomb interaction, and 
H e -ph=gQ p+P 2 /2+V (Q) . Here g is the electron-phonon cou- 
pling constant, p=n^+ni, Q is normal coordinate of the os- 
cillation of a caged atom, P is the corresponding canonical 
momentum, and V(Q) indicates the cage potential. Note that 
the reduced mass of the vibration is set as unity. The poten- 
tial is given by V(Q)=u> 2 Q 2 /2+fc4<2 4 +fc 6 Q 6 , where u is the 
phonon energy and we consider fourth- and sixth-order anhar- 
monicity in the potential. Using the phonon operator a defined 
through Q=(a + a^)/\/2w, we obtain 

ifeph = ujy^a(a + a')p + uj(a)a + 1/2) 

+ (3 4 uj(a + a f ) 4 + 6 w(a + a f ) 6 , 



(2) 



where the non-dimensional electron-phonon coupling con- 
stant a is given by a=g 2 /(2u 3 ) and non-dimensional an- 
harmonicity parameters /Jj and [3q are, respectively, given 
by /?4=fc4/(4w 3 ) and (3q=}iq / (%<ju a ) . By using these non- 
dimensional parameters, we express the potential as 

auj(q 2 + 16a/3 4 g 4 + 64a 2 /3 6 g 6 ), 



V(q) 



(3) 



where q is the non-dimensional length, given by q=Quj 2 /g. 

In this paper, we fix f3e as (3g=10~ 5 and change /3 4 in the 
negative side in order to examine the effect of potential shape. 
Note that for /?4>0, there is no drastic change in the potential 
shape. As for a, we fix a=2 to consider the strong electron- 
phonon coupling. In the calculations, we define the phonon 
basis as \£)=(af) e \0) /VH, where I is the phonon number and 
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Fig. I. (a) Change of the potential shape due to /?4 for «=2 and /36=10 _ . 
Broken curve denotes the quadratic potential, (b) Average value of q vs. 
P4 for p=l and 2. (c) Phonon-mediated attraction Up^/w vs. P4. 



|0) is the vacuum state. The phonon basis is truncated at a 
finite number, which is set as 2000. 

In Fig. 1(a), we depict V(q)/(au>) in order to see how 
the potential shape is changed when 04 is decreased in the 
negative side. For /?4=0, there is a single minimum at q=Q. 
In comparison with the case of quadratic potential, the an- 
harmonic potential becomes steep for large q, while there 
is no significant difference among the potentials in the re- 
gion near q=0. Thus, the potential in this parameter region is 
called the on-center type. When (3^ is decreased, shoulder-like 
structure gradually appears and the potential bottom becomes 
flat and wide. In the present definitions of parameters, poten- 
tial minima at |g|^0 appear for /3 4 <-\/3/3 6 /4=-0.00274. 
This shape of the potential is called the rattling type, since 
the caged atom is expected to oscillate with large ampli- 
tude in spite of the structure with shallow minima in the 
wide bottom of the potential. For further decrease of /3 4 with 
pi<—\/P&=— 0.00316, the minimum values of V(q) at q^Q 
are smaller than that at the origin. The shape of the potential 
in this region is called the off-center type in the sense that 
the atom oscillates around at off-center positions due to deep 
potential minima at |g|^0. 

Here we summarize the three types of the potential shape: 
(i) The on-center type for - t/3/? 6 /4</?4<0, (ii) the rattling 
type for -VS<ft<-\/W^' ar, d ( m ) the off-center type 
for (3i<—\f^Q. When we include the electron-phonon cou- 
pling, the boundary values will be changed, but the above 
three types still characterize the potential shape. 

In order to see the effect of the potential shape on the elec- 
tron state, let us evaluate the average displacement (q) for p=l 
and 2. As shown in Fig. 1(b), for the case of p=l, when we 
decrease 04 from zero, first (q) is gradually decreased, since 
the bottom of the on-center type potential is not so widened 
by 04. Around at /3 4 «— 0.0025, corresponding to the region 



of the rattling type, (q) is rapidly decreased. When we further 
decrease 04, (q) is gradually decreased for (5a<— 0.003 in ac- 
cordance with the shift of potential minima. For p=2, (q) is 
smoothly decreased in the whole region of /34<0. Note that in 
the regions for on-center and rattling types, we observe signif- 
icant difference between the results for p=l and 2. Since the 
electron-phonon coupling for p=2 is virtually stronger than 
that for p=l, the positions of potential minima for p=2 are 
larger than that for p=l. However, in the region of the off- 
center type, the potential is heavily deformed by large nega- 
tive value of (3 4 even without the electron-phonon coupling. 
Thus, the position of the potential minima is insensitive to the 
electron number for /?4< — 0.003. 

The difference in the dependence of the change of the po- 
tential shape on electron number can be significantly found 
in the effective attraction between electrons mediated by 
phonons. When we define the attraction as — f7 p h, the mag- 
nitude of the attraction U p h is given by 

U ph = 2E^-(E^+E^), 



(4) 



where is the ground-state energy of H ep h for the electron 
number p. We note that the effect of U is not included in E p . 
Note also that U p \ l =2auj at half-filling for harmonic phonon. 

In Fig. 1(c), we plot Up^/us as a function of 04. At /9 4 =0, 
Uph/tu is about 3.7, which is close to 2a for a=2. When we 
decrease 04, <7 p h is increased, since the degree of the poten- 
tial deformation for p=2 is relatively larger than that for p=l. 
When the bottom of the potential becomes flat and wide, the 
phonon density of states at low energies are increased. Then, 
the polaronic effect is enhanced, leading to the increase of the 
attraction. However, in the region of the off-center type, there 
is no significant difference in the degree of the potential defor- 
mation between the cases for p=l and 2, as understood from 
Fig. 1(b). Then, <7 p h is decreased when we further decrease 
04 in the region of the off-center type potential. 17 

From the above result, we can understand that the Coulomb 
interaction is effectively reduced for the potential shape of 
the rattling type. For a reasonable choice of the parameter as 
U /oj=\Q, we actually find a region with effective attraction 
around at /3 4 =— 0.0025. This suggests that the Kondo temper- 
ature should be enhanced in such a region. In the following, 
let us focus on the Kondo effect for -0.003</3 4 <-0.002. 

Now we include the hybridization between localized and 
conduction electrons. Since we consider localized oscillation 
of the caged atom, we ignore the coupling between phonons 
and conduction electrons, Then, the Hamiltonian is given by 

H = ^2£kc{ a c ka + ^2(Vcl a d a + h.c.) + e d p+H loc , (5) 

ha ka 

where e k denotes the dispersion of conduction electron, c k a is 
an annihilation operator of conduction electron with momen- 
tum k and spin a, V is the hybridization between conduc- 
tion and localized electrons, and Ed indicates a local electron 
level. In the following, we always adjust Ed to consider the 
half-filling case, i.e., (p)=l, even though we do not mention 
explicitly. The energy unit is chosen as a half of the conduc- 
tion bandwidth, which is expressed by D. Hereafter we fix 
V/D=0.25 and U/D=2. Typically, D is in the order of 1 eV. 

In order to clarify electronic properties of H, we evaluate 
entropy Si mp , specific heat Cj mp , and susceptibility by using a 
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Fig. 2. (Color online) (a) T\s and T\ c vs. temperature, (b) Simp 
Ci mp vs. temperature, (c) Kondo temperatures vs. /?4. 
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Fig. 3. (Color online) (a) Sommerfeld constant 7 vs. (3^ for several values 
of u). (b) 7 in the unit of mJ/K 2 -mol vs. magnetic field for D=5000 K. 



numerical renormalization group (NRG) method. 18 The loga- 
rithmic discretization of the momentum space is characterized 
by a parameter A and we keep M low-energy states for each 
renormalization step. Here we set A=2.5 and Af=5000. Spin 
and charge susceptibilities are, respectively, given by 

1 ^ e -E n /T _ e -E m /T 

Xs = - > J £ 5 \(n\a z \m)\ , 



Z ^ 

n,m 

1 x -> e 



-E n /T 



e -E m /T 



(6) 



E,, 



E n 



\(n\(p-( P ))\m)f 



where E n is the eigenenergy for the n-th eigenstate | n) of H , 
T is a logarithmic temperature defined as T=DA~( N ~ 1 * > / 2 , 
Z is the partition function given by Z=^2 n e~ E '^ T , {p)= 
E„ e~ E -/ T (n\p\n), and a z =n T -n x . 
In Figs. 2(a) and 2(b), we show typical results for sus- 
ceptibility, entropy, and specific heat for /?4=— 0.0022. Note 
that this is the region with weak effective attraction, given 
by U— C/ p h=— 0.075D. Thus, the spin susceptibility is rapidly 
suppressed due to the effect of on-site attraction, while the 
charge susceptibility is increased without the slight decrease 
near T / D=\. When the spin susceptibility is suppressed, spin 
entropy log 2 is released, leading to a peak in the specific heat. 

(s) 

This peak temperature is defined as T-^ ' . 

After the release of the spin entropy, we observe a narrow 
plateau of log 2, corresponding to the degenerate region char- 
acterized by vacant and double occupied states. This entropy 
of log 2 is gradually released and it eventually goes to zero 
around at T/ D^10~ 3 , where we can see another peak in the 

(c) 

specific heat. We define this temperature as . 

Changing the values of /?4, we repeat the calculations and 
plot the Kondo temperatures in Fig. 2(c). As understood from 
Fig. 2(b), we can observe a couple of peaks in the specific 
heat, corresponding to the vanishment of spin and charge sus- 
ceptibilities. Thus, we carefully check the behavior of Txs 
and Txc, in order to specify the origin of the entropy release, 



leading to the two kinds of the Kondo temperature. 

(si 

In the region of the on-center type potential, is not so 
affected by @4 and we cannot find a high-temperature peak 

(c) 

corresponding to . However, when we approach the re- 

(s) (c) 

gion of the rattling-type potential, is increased and 
appears due to the reduction of the Coulomb interaction. 
Eventually, the effective interaction turns out to be attractive 

(c) (s) 

and T-^ ' becomes lower than ' . When we further decrease 
/?4, we enter the region of the off-center type and the phonon- 
mediated attraction is rapidly decreased, as found in Fig. 1(c). 
The peak corresponding to the entropy release of charge de- 
gree of freedom is changed to just a shoulder-like structure, 
and it finally disappears. 

In general, when the Kondo temperature is increased, we 
can expect that 7 is magnetically robust. Roughly speaking, 
the electronic properties will not be significantly affected by a 
magnetic field H smaller than fceTk/ Mb> where p^ is a Bohr 
magneton. Thus, in the region of the rattling-type potential 
with enhanced Tk, we can expect that 7 is magnetically ro- 
bust. On the other hand, in the Kondo problem, it has been 
known that 7 is in proportion to 1/Tk- 16 As we will describe 
later, the behavior of 7 as a function of f3i is well explained 
by that of 1/Tk. When the Coulomb repulsion is weakened 
by the phonon-mediated attraction, the Kondo temperature is 
increased and thus, 7 is decreased. Namely, the magnitude of 
7 is not expected to be large in the region of the rattling type. 

However, in the simple discussion based on the relation of 
7^1/Tk, it is unclear whether the magnitude of 7 is large or 
not. Thus, we actually evaluate 7, given by 7=Ci mp /T at low 
temperatures. Here 7 is estimated at the lowest temperature at 
which we can reach in the present NRG calculation. 

In Fig. 3(a), we show 7 vs. (3^ for several values of uj. 
For u> j D=0.01, in which the adiabatic approximation is con- 
sidered to be valid, the value of 7 is not so sensitive to /?4, 
even if the potential shape is changed by ^4. In particular, for 
/?4>— 0.002, the values of 7 are similar to that of the Ander- 
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son model without phonons. The deviation can be found only 
in the region of the potential of the off-center type. 

When we increase the value of u, the effect of anti- 
adiabaticity appears. For w=0.2, 7 decreases significantly 
for — 0.003</?4<— 0.002, which is, roughly speaking, corre- 
sponding to the region of the rattling-type potential. Note that 
this behavior of 7 is well reproduced by 1 / . When we fur- 
ther increase u>, the region of the suppressed 7 becomes wide 
and for uj/D=0.5, such a region appears for /3 4 >— 0.003. 
Note that very large values of 7 can be found only in the re- 
gion of the off-center type potential. 

Now we turn our attention to the effect of a magnetic field. 
In Fig. 3(b), we show 7 vs. magnetic field H along the z-axis 
for -0.003</? 4 <-0.002 with a=2, u/D=Q.% and /3 6 =10~ 5 . 
In this calculation, we add the Zeeman term H mag =ii-Qa z H 
to the Hamiltonian (5), where the g-factor is set as 2. To eval- 
uate the values of 7 in the unit of mJ/K 2 -mol, we set Z?=5000 
K. For /3 4 deviated from the rattling-type potential, the mag- 
nitude of 7 at H=0 is enhanced, but it decreases with the 
increase of H, as found in the curves for /?4=— 0.002 and 
—0.003. On the other hand, when we take /J4 corresponding to 
the rattling-type potential, 7 is actually magnetically robust, 
but the magnitude of 7 is not so large. However, when we 
consider the parameter region just between the on-center and 
the rattling types (/?4=— 0.0021) or between the off-center and 
the rattling types (/?4=— 0.0028), it seems to be possible to re- 
produce magnetically robust 7 in the order of 100 mJ/K 2 -mol. 
These regions are corresponding to the cases at which the 
Coulomb repulsion is canceled by phonon-mediated attrac- 
tion. Thus, on the basis of a concept of the Kondo effect in the 
electron-phonon system, we propose that magnetically robust 
heavy fermion behavior is understood by the competition be- 

(s) (c) 

tween spin and charge fluctuations at the region of '. 

Here we provide a couple of comments, (i) It has been 
pointed out that the competition between the Coulomb in- 
teraction and phonon-mediated attraction plays an important 
role in strongly coupled electron-phonon system. 19 24 Here 
we control such a competition by anharmonicity to explain 
magnetically robust heavy fermion phenomenon, (ii) We have 
considered the one-dimensional potential for simplicity, but 
in actuality, the motion of the caged atom occurs in a three- 
dimensional potential. In particular, there exists a mode in 
which the caged atom moves around potential minima. 1 Con- 
cerning this issue, the present author has analyzed the An- 
derson model coupled with Jahn-Teller phonons. 9 Then, we 
have proposed an intriguing phenomenon of chiral Kondo ef- 
fect due to the rotational degree of freedom of Jahn-Teller 
phonon. 10 On the basis of this effect, the magnetically ro- 
bust heavy fermion phenomenon may be explained for the 
off-center type potential. It is one of future issues. 

Finally, we briefly discuss the effect of orbital degree of 
freedom, which has been ignored in this paper. Although the 
motivation of this work is magnetically robust heavy fermion 
phenomenon observed in Sm-based filled skutterudites, we 
have not considered / electrons in the model. In /-electron 
systems with orbital degeneracy, multipole degree of free- 
dom becomes active due to strong spin-orbit coupling. In this 
sense, in order to discuss electronic properties of actual Sm- 
based filled skutterudites, it is important to consider multi- 
pole degrees of freedom. 25 The multi-orbital Anderson model 



coupled with Jahn-Teller phonons has been also analyzed, 
but further works are required along with this research di- 
rection, in particular, in order to understand the difference in 
SmOs4Sbi2 and SmFe4Pi2. It is also expected that more sta- 
ble magnetically robust heavy fermion state is realized when 
orbital fluctuations are also comparable with spin and charge 
fluctuations. It is another future problem. 

In summary, we have discussed the Kondo effect of the 
Anderson model coupled with local anharmonic phonon. In 
order to understand the magnetically robust heavy fermion 
phenomenon, we have emphasized a key issue of the shape 
of the cage potential. In particular, we have found that 7 is 
both large and magnetically robust, when the potential shape 
is deformed from the rattling type, in which spin and charge 
fluctuations are comparable to each other. 
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